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Abstract. When applied to Be/X-ray binaries, the viscous decretion disc model, which can successfully account 
for most properties of Be stars, naturally predicts the truncation of the circumstellar disc. The distance at which 
the circumstellar disc is truncated depends mainly on the orbital parameters and the viscosity. In systems with 
low eccentricity, the disc is expected to be truncated at the 3:1 resonance radius, for which the gap between the 
disc outer radius and the critical lobe radius of the Be star is so wide that, under normal conditions, the neutron 
star cannot accrete enough gas at periastron passage to show periodic X-ray outbursts (Type I outbursts). These 
systems will display only occasional giant X-ray outbursts (Type II outbursts). On the other hand, in systems 
with high orbital eccentricity, the disc truncation occurs at a much higher resonance radius, which is very close 
to or slightly beyond the critical lobe radius at periastron unless the viscosity is very low. In these systems, disc 
truncation cannot be efficient, allowing the neutron star to capture gas from the disc at every periastron passage 
and display Type I outbursts regularly. In contrast to the rather robust results for systems with low eccentricity 
and high eccentricity, the result for systems with moderate eccentricity depends on rather subtle details. Systems 
in which the disc is truncated in the vicinity of the critical lobe will regularly display Type I outbursts, whereas 
those with the disc significantly smaller than the critical lobe will show only Type II outbursts under normal 
conditions and temporary Type I outbursts when the disc is strongly disturbed. In Be/X-ray binaries, material 
will be accreted via the first Lagrangian point with low velocities relative to the neutron star and carrying high 
angular momentum. This may result in the temporary formation of accretion discs during Type I outbursts, 
something that seems to be confirmed by observations. 
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1. Introduction 

Be/X-ray binaries are X-ray sources composed of a Be star 
and a neutron star. The high-energy radiation is believed 
to arise owing to accretion of material associated with the 
Be star by the co mpact object (see Negueruela 1998; see 
also Bildsten et al. |l997|) . 



A "Be star" is an early- type non-supergiant star, which 
at some time has shown emission in the Balmer series lines 
(Slettebak 1988, for a review). Both the emission lines 
and the characteristic strong infrared excess when com- 
pared to normal stars of the same spectral types are at- 
tributed to the presence of circumstellar material in a disc- 
like geometry. The causes that give rise to the disc are not 
well understood. Different mechanisms (fast rotation, non- 
radial pulsation, magnetic loops) have been proposed, but 
it is still unclear whether any of them can explain the ob- 
served phenomenology on its own. The discs are rotation- 



ally dominated and motion seems to be quasi-Keplerian 
(Hanuschik 1996). However, some kind of global outflow 
is needed to explain the X-ray emission from Be/X-ray 
binaries (Waters et al. 198^ ). 

Some Be/X-ray binaries are persistent X-ray sources 
(see Reig & Roche 1999D , displaying low luminosity {L^ 
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lO'^^erg s^^) at a relatively constant level (varying by up to 
a factor of ~ 10). On the other hand, most known Be/X- 
ray binaries (though this is probably a selection effect) 
undergo periods in which the X-ray luminosity suddenly 
increases by a factor ^ 10 and are termed Be/X-ray tran- 
sients. 

Be/X-ray transients fall within a relatively narrow area 
in the Por b/^sp in diagram (see Corbet 1986; Waters & van 
Kerkwijk |l989| , indicating that some mechanism must be 
responsible for the correlation. Those systems with fast- 
spinning neutron stars do not show pulsed X-ray emission 
during quiescence (though non-pulsed radiation could be 
caused by accretion on to the magnetosphere) because of 
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the centrifugal inhibition of accretion (Stella et al. 19861 ). 
Systems with more slowly rotating pulsars show X-ray 
emission at a level ^ 10"^^ erg s^^ when in quiescence. 
Transients show two different kinds of outbursts: 

— X-ray outbursts of moderate intensity (Lx ~ 10 
10'^^ erg s~^) occurring in series separated by the or- 
bital period (Type I or normal), generally (but not 
always) close to the time of periastron passage of the 
neutron star. In most cases, the duration of these out- 
bursts seems to be related to the orbital period. 

— Giant (or Type II) X-ray outbursts {L^ ^ 10^^ erg s~^) 
lasting for several weeks or even months. Generally 
Type II outbursts start shortly after periastron pas- 
sage, but do not show any other correlation with or- 
bital parameters (Finger & Prince 1997). In systems 
like 4U 0115-1-63 the duration of the Type II outbursts 
seems to be to some degree correlated with its peak 
intensity, but in A 0535-1-262 Type I outbursts may be 
as long as much brighter Type II outbursts (Finger et 
al. p96a|). 



2. Radial outflows vs. quasi-Keplerian discs 

Attempts at modelling the X-ray luminosities of Be/X- 
ray transients during outbursts have made use of a simple 
wind accretion model, in which the neutron star accretes 
from a relatively fast radial outflow. The disc of the central 
Be star is supposed to have a power-law density distribu- 
tion 



(1) 



P{r) = PO ( ^ 



where po is the density at the stellar surface and i?* is the 
radius of the star. This results in a power velocity law of 
the form 

n-2 



v{r) = vo 



(2) 



where the values of vq and n have to be determined obser- 
vationally (see Waters et al. 1989 and references therein). 
The rotational velocity of the outflow takes the form 



(3) 



with 0.5 < a < 1 (respectively the Keplerian case and 
conservation of angular momentum). 

Accretion is considered to follow the classical Boiidi- 
Hoyle-Littleton (BHL) approximation. The most impor- 
tant parameter, the relative velocity between the outflow 
and the neutron star, can be written as 



{V - Vrad) + (Wrot ~ Kot)' 



(4) 



where K-ad and K-ot are the radial and tangential compo- 
nents of the orbital velocity of the neutron star. Material 
is supposed to be accreted when it is within a capture 
radius defined as 



where Mx is the mass of the neutron star and G is the 
gravitational constant. The X-ray luminosity in the BHL 
approximation can then be expressed as 

(6) 



Lx = AttG^M^K 



-4 T-, -3 



where i?x is the radius of the neutron star and Fm — pvrd 
is the mass flow. In order to explain the wide range of 
observed X-ray luminosities, large changes in the value 
of the radial velocity have to be invoked. For example. 



Waters et al. (1989) deduced that the relative velocity was 
Wj-oi w 300 kms^^ during a Type I outburst of V 0332+53 
in 1983, while it was ^ 100 km s^^ during a Type II out- 
burst in 1973. 

The use of the BHL approximation implies a number 
of simplifying assumptions which are not always easy to 
justify. For example, it neglects any effect of the mass- 
losing star, which for periastron distances of ~ 10 i?* and 
mass ratios q ^ 0.1 seems to be an excessive simplification. 
In addition, while the use of the accretion radius formalism 
is adequate for an accreting object immersed in a medium, 
its application to an outflow characterised by a relatively 
small scale-height (H), such as a Be disc, is dubious. 

Moreover, when the low outflow velocities required to 
explain Type II outbursts are considered, the formalism 
breaks down completely, since the capture radius becomes 
far too large to have any physical meaning. For example, 
for Wrci ^ 20 kms-\ ~ 9.6 x 10" m « 1400 Rq, which 
is one order of magnitude larger than the binary separa- 
tion. A crude way round this problem is to consider that 
the radius of the effective Roche lobe of the neutron star 

should be used instead of rc whenever the calculated 
value is larger than (e.g., Ikhnasov 2001). In spite of 



(5) 



these shortcomings, the model has been repeatedly used 
in an attempt to mod el ligh tcurves of B e/X-r ay binaries 
(Raguzova & Lipunov 1998 Reig et al. 1998| ) with only 
moderate success. 

Beyond the purely formal aspects, one obvious diffi- 
culty for the model is the fact that many Be/X-ray bina- 
ries (e.g., A 0535-1-262) show low-luminosity X-ray emis- 
sion when they are not in outburst. It is believed that 
all Be/X-ray binaries for which centrifugal inhibition of 
acc retion is not effective display this emission. Motch et 
al. (1991) detected A 0535-1-262 on several occasions at lu- 
minosities of 2 X 10^^ ergs~^. In order to explain this 
luminosity within the framework described above and tak- 
ing into account that optical and infrared observations do 
not show any sign of the large variations that would be 
associated with a change of several orders of magnitude 
in the density of material, enormous relative velocities (of 
the order of ~ 10'* kms~^) are needed. 

One further complication comes from the fact that 
Be/X-ray binaries spend most of their time in the qui- 
escent state described in the previous paragraph and only 
occasionally show series of outbursts. The model does not 
offer any explanation as to why there could be a change 
from quiescence to outburst, unless again very large and 
sudden changes in the density and velocity of the flow are 
assumed. Given the changes in relative velocities needed 
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to account for the observed range of X-ray luminosities 
and the lack of any physical mechanism that could ex- 
plain them, it is clear that direct accretion from a wind- 
like outflow is not the best approximation to the way in 
which matter is fed on to the neutron star. 

But the major objection to the model is simply the 
fact that there is no observational evidence whatsoever 
supporting the existence of such fast outflows. All obser- 
vations of Be stars imply bulk outflow velocities smaller 
than a few kms~^ (Hanuschik 2000). The evidence for 



rotationally dominated quasi-Keplerian di scs a round Be 
stars is overw helm ing (see Hanus chik et al. 1996| ; Hummel 
& Hanuschik 1997 ; Okazaki 1997 ), especially owing to the 
success of the one-armed global oscillation model to ex- 
plain V/R variabilit y in the em ission lines of Be stars 
(Kato |1983| ; Okazaki |l99l| |l996| ; Papaloizou et al. |l992 ; 
Hummel & Hanuschik 1997 ). 

Therefore it seems necessary to attempt an explana- 
tion of the outburst behaviour of Be/X-ray binaries that 
does not imply large outflow velocities. 

3. The viscous disc model 

Whatever the mechanism originating the Be phenomenon, 
the model which at present appears more applicable to 
explaining the discs surrounding Be stars is the viscous 
decretion disc model (Lee et al. 1991; see Porter 1999 and 
Okazaki 2001 for detailed discussion). In this scenario, an- 



gular momentum is transferred from the central star by 
some mechanism still to be determined (perhaps associ- 
ated with non-radial pulsations) to the inner edge of the 
disc, increasing its angular velocity to Keplerian. Viscosity 
then, operating in a way opposite to an accretion disc, 
conducts material outwards. In this scenario, material in 
the disc moves in quasi-Keplerian orbits and the radial 
velocity component is highly subsonic until the material 
reaches a dista nce m uch larger than the line-emitting re- 
gions (Okazaki 2001). The outflow is very subsonic for the 
distances at which neutron stars orbit in close Be/X-ray 
transients (v{r) < 1km s^"'^) and still subsonic for the or- 
bital sizes of all Be/X-ray binaries for which there is an 
orbital solution. The viscous decretion disc model success- 
fully accounts for most of the observational characteristics 
of Be discs. 



Negueruela & Okazaki (2001, henceforth Paper I) have 
modelled the disc surrounding the Be primary in the 
Bc/X-ray transient 4U 0115+63 as a viscous decretion 
disc and found that the tidal interaction of the neutron 
star naturally produces the truncation of the circumstel- 
lar disc, as it does for accretion discs in close binaries 



(Paczyhski 1977). The result of Negueruela & Okazaki 
(2001) is in agreement with the results of Reig et al. 
(1997), who showed that there is a correlation between 
the orbital size and the maximum equivalent width of Ha 
ever observed in a system. Even though it is clear that the 
equivalent width of Ha is not an effective measurement of 
the size of the disc owin g to several effects (see, for ex- 
ample, Negueruela et al. 1998 ), the maximum equivalent 



width ever observed becomes a significant indicator if the 
system has been monitored during a period which is long 
in comparison with the typical time-scale for changes in 
the disc (which, if viscosity is dominant, should be only 
a few months). Therefore the result of Reig et al. (1997) 
clearly indicates that the neutron star has some sort of 
effect on the size of the disc. 

In this paper we apply the model presented in Paper I 
to several Be/X-ray transients for which orbital solutions 
exist and investigate how the truncation radius depends 
on different orbital parameters. 

4. Model description and limitations 

The model developed in Paper I describes a binary system 
in which a primary Be star of mass M^, and radius i?* is 
orbited by a neutron star of mass = 1.4 Mq which 
moves in an orbit of eccentricity e and period Porb- The 
Be star is assumed to be surrounded by a near-Keplerian 
disc which is primarily governed by pressure and viscosity. 
For simplicity, the disc is assumed to be isothermal and 
Shakura-Sunyaev's viscosity prescription is adopted. 

In such a disc, angular momentum is added to the 
disc by the viscous torque, whereas it is removed from the 
disc by the resonant torque exerted by the neutron star 
companion, which becomes non-zero only at radii where 
the ratio between the angular frequency of disc rotation 
and the angular frequency of the mean binary motion is 
a rational number. As a result, the disc decretes outward 
owing to the transfer of angular momentum by viscosity 
until the resonant torque becomes larger than the viscous 
torque at a resonant radius. 

Therefore, the criterion for the disc truncation at a 
given resonance radius is written as 



Tvis + Tres < 0, 



(7) 



where Tvis and Tj-cs are the viscous torque and the resonant 
torque, respectively. 

The viscous torque Tvis is written as 



STraGAf* err 



H 



(8) 



(Lin & Papaloizou 1986 ), where a is the surface density 
of the disc, a is the Shakura-Sunyaev viscosity parameter 
and H is the vertical scale-height of the disc given by 



H 



Vk{R*) \R 



1/2 



(9) 



for the isothermal disc. Here, Cs is the isothermal sound 
speed and Vk{R*) is the Keplerian velocity at the stellar 
surface. In the systems we will discuss later, Cs/Vk{R*) 
ranges 3.1 - 4.1 • IQ-^{Ta/T^sY/^ , where Td and T^s are 
the disc temperature and the effective temperature of the 
Be star, respectively. Note that the viscous torque is pro- 
portional to the disc temperature. 

The resonant torque Ti-cs is calculated by using 
Goldreich & Tremaine's ( 1979 , 1980 ) torque formula, after 
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decomposing the binary potential $ into a double Fourier 
series as 



$(r,6',z) 



[r^+rl-2rr2 cos(e' - /)]i/2 



H —cos{e - / 

(/-m/ exp [1(7716* - l^Bt)], 



(10) 



where r2 is the distance of the neutron star from the pri- 
mary, / is the true anomaly of the neutron star, m and 
/ are the azimuthal and time-harmonic numbers, respec- 
tively, and Qb — [G(M, + M^)/a'^]^/^ is the mean motion 
of the binary with semimajor axis a. The pattern speed of 
each potential component is given by f2p — {l/m)QB- The 
third term in the right hand side of the first equation is the 
indirect potential arising because the coordinate origin is 
at the primary. 

For each potential component, there can be three kinds 
of resonances, i.e., the outer and inner Lindblad reso- 
nances at radii [{m ± l)//]^/'^(l -I- q)~^^'^a, where flp — 
± K/m, and a corotation resonance (CR) at the ra- 
dius (m/Z)^/^(l -I- g)~^/^a, where fip = il. Here, k is 
the epicyclic frequency and the upper and lower signs 
correspond to the outer Lindblad resonance (OLR) and 
inner Lindblad resonance (ILR), respectively. In circum- 
stellar discs, however, the resonant torque from the inner 
Lindblad resonance, which is given by 



{Trnl) 



m/JILR, 



m(m — l)7r^a-(A -|- 2m)^0: 



ral 



(11) 



where A = d In (pmi /d In r, always dominates the reso- 
nant torques from the corotation resonance and the outer 
Lindblad resonance. The resonant torque at a given reso- 
nance radius is then given by 

m' I m ' l' 

(12) 



ml 



ILR- 



Since high-order potential components contribute little to 
the total torque, the summation in Eq. ( p^ is safely taken 
over several lowest-order potential components which give 
the same radius. 

For a given set of stellar and orbital parameters and 
the disc temperature, criterion (^ at a given resonance is 
met for a smaller than a critical value Ofcrit and we assume 
that the disc is truncated at the resonance if a < acrit- 

In Paper I it was shown that for any subsonic outflow, 
the drift time-scale Tdrift Ar/vr M~^{Ar/H)fl~^ 
was considerably longer than a typical truncation time- 
scale Ttrunc ~ (a/Q!crit)Tvis ~ a~]^{Ar / H)'^n~^ , whcrc Ar 
is the gap size between the truncation radius and the ra- 
dius where the gravity by the neutron star begins to dom- 
inate, and Vr and A4r are the radial velocity and Mach 



number, respectively. As a consequence, the tidal and res- 
onant interaction with the neutron star led to disc trun- 
cation. 

One important consequence of the above is that the 
discs surrounding the primaries in Be/X-ray binaries can- 
not reach a steady state. Most of the outflowing mate- 
rial loses angular momentum and falls back towards the 
central star. As a consequence of the interaction with the 
material which is coming outwards from the inner regions, 
it is likely that the disc becomes denser and the density 
distribution in the radial direction becomes flatter with 
increasing time. We note that we do not expect the trun- 
cation effect to be one hundred per cent efficient. This 
is not only owing to theoretical considerations (see be- 
low), but probably required by the existence of pulsed 
low-luminosity X-ray emission during quiescence. 

In the formulation above only torques integrated over 
the whole orbit are considered. Given the large eccentric- 
ities observed in Be/X-ray binaries (mostly larger than 
0.3 and sometimes approaching 0.9), the gravitational ef- 
fect of the neutron star is very strongly dependent on the 
orbital phase. This means that in such systems the disc ra- 
dius and the truncation radius are also phase-dependent. 
The disc would shrink at periastron, at which the trunca- 
tion radius becomes smallest, while it would spread when 
the neutron star is far away and its gravitational effect is 
not felt so strongly. The spread of the disc would continue 
until the truncation radius becomes smaller than the disc 
radius at some phase before the next periastron passage. 
This variation in disc radius would be larger for a larger 
eccentricity. For a longer orbital period, the effect would 
be yet stronger. Given that, for some systems considered 
(see Section |^), the model truncation radius is close to 
the critical lobe radius at periastron, this will provide a 
mechanism by which disc material can reach the neutron 
star (see also Fig. ||for scenarios for two families of Type I 
outbursts). 

5. Modelling the Be/X-ray binaries 

Table |l] shows a list of known X-ray binaries which have 
exhibited outbursting behaviour or some sort of orbital 
modulation in their X-ray lightcurve. The top panel con- 
tains Be/X-ray binaries with an identified Be counterpart 
which have displayed Type I X-ray outbursts, i.e., a series 
of outbursts separated by their orbital period (note that 
in many cases there is no orbital solution for the system 
and the recurrence period of the outbursts is taken to be 
the orbital period). The middle panel shows other sys- 
tems without identified optical counterparts whose X-ray 
behaviour marks them out as Be/X-ray binaries. Finally 
the bottom panel contains a few other Be/X-ray binaries 
with identified counterparts whose X-ray behaviour devi- 
ates slightly from what is considered typical (and which 
will be discussed individually in Section 

In this section, we apply our model to systems for 
which exact orbital solutions have been deduced from the 
analysis of Doppler shifts in the arrival times of X-ray pho- 
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tons. They include five Be/X-ray binaries in the top panel 
(V 0332+53, A0535+262, GRO J1008~57, 2S1417-624, 
and EXO 2030+375) and a likely Be/X-ray binary in the 
middle panel (2S 1845—024). These systems are discussed 
individually in the following subsections. The case of the 
Be/X-ray transient 4U 0115+63 has been carefully dis- 
cussed in Paper I and therefore it will not be included 
here. 

As mentioned in the previous section, we adopt 
Shakura-Sunyaev's viscosity prescription, in which the vis- 
cosity parameter a is a free parameter, and assume the 
Be disc to be isothermal. In what follows, we adopt Td = 
^Teff. Note that our assumption of the disc temperature 
is consistent with the results by Millar and Marlborough 
( |1998| , |1999D , who computed the distribution of the disc 
temperature within 100 i?* around the BO star 7 Gas and 
the B8-9 star 1 Del by balancing at each position the 
rates of energy gain and energy loss and found that the 
disc is roughly isothermal at a temperature about half the 
effective temperature of the star. 

When adopting a particular model for a system, the 
main source of uncertainty comes from the choice of mass 
for the primary, a parameter which can only be guessed 
from the spectral type. We find two main difficulties. In 
some cases, the spectral type of the primary is not well 
determined. Moreover, there is some evidence that fast ro- 
tators may be moder ately over-luminous for their masses 
(see Gies et al. 1998 ). For this reason, in general we have 
taken masses slightl y low er than those given in the calibra- 
tion of Vacca et al. ( 1996| ). The spectral distribution of pri- 
maries of Be/X-ray binaries (Negueruela 1998 ) is strongly 
peaked at BO. Therefore for systems without exact deter- 
mination of the spectral type, we have calculated models 
corresponding to BOV and BOIII primaries. In any case, 
our results show that the exact mass of the primary is 
not one of the main factors in the X-ray behaviour of the 
sources. 

5.1. V 0332+53 

This transient pulsar has an orbital period Porb = 34.25 d 
with a relatively low eccentricity of e = 0.31 (Stella et 
al. 1985). The optical component of this system is an un- 
evolved star in the 08-9 range (Negueruela et al. 19991 ). 
For 08.5V stars Vacca et al. ( |1996D give an spectroscopic 
mass M, — 23.6 Mq and a theoretical mass Af, — 28 Mq. 
Here we will assume as a conservative model the lower 
Hmit of M, = 20 Mq and i?* = 8.8 R(^. The mass func- 



tion f{M) — 0.1 then implies sini — 0.17, resulting in an 
orbital separation a ~ 130 i?0 (rpcr ~ 10 i?,). 

This system has been rarely detected in X-rays. A 
Type II out burst was observed in 1973 (see Negueruela 
et al. 1999 for references). Ten years later, it was ob- 
served during a series of three Type I outbursts. Finally 
it was observed during another Type II outburst in 1989. 
Between 1991 and 2000 it has not been detected by either 
the BATSE instrument on board the ComptonGRO satel- 
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Fig. 1. Gritical values of a at some resonance radii for 
systems discussed in the text. Annotated in the figure are 
the locations of the n : 1 commensurabilities of disc and 
mean binary orbital frequencies. Td = ^Tefi is adopted for 
all models. For other disc temperatures, acrit should be 
multiplied by a factor of Tcii/2T^. 

lite or the All Sky Monitor on board RossiXTE and it is 



believed to be in a dormant state (Negueruela et al. 1999). 

In Fig. 0, we plot acrit at the n : 1 resonance radii for 
those Be/X-ray binaries which will be discussed in this sec- 
tion. The resonant torques at the n : 1 radii are stronger 
than those at radii with other period commensurabilities 
located nearby. Note that we have adopted a particular 
disc temperature, Td = ^Tcff, for each stellar model. For 
other disc temperatures, acrit should be multiplied by a 
factor of Tcff/2Td, taking account of the fact that the vis- 
cous torque is proportional to Td. 

Fig. shows that the Be disc in V 0332+53 is expected 
to be truncated at the 3:1 resonance radius {rt/a ~ 0.47, 
where rt is the truncation radius) for 0.099 ^ a ^ 0.60 
and at the 4:1 resonance radius (n/a ~ 0.39) for 0.019 ^ 
a <, 0.099. 

It is interesting to see how close the truncation radius is 
to the size of the critical lobe at periastron. Fig. || shows 
orbital models for the systems discussed in this section. 
The potential -0 describing the effects of the gravitational 
and centrifugal forces on the motion of test particles or- 
biting the Be star is given by 



V'(r, 61, z) = $(r, e, z) - -VL^{r)r'^, 



(13) 



where $ is the potential defined by Eq. (|l^). Also shown is 
the distance scale corresponding to O.lCgPorb- This scale 
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Table 1. List of known X-ray binaries to which our model could be applied. The top panel contains Be/X-ray binaries 
with an identified counterpart which have displayed Type I X-ray outbursts. The middle panel shows other systems 
without identified optical counterparts likely to be Be/X-ray binaries. The bottom panel contains some other Be/X- 
ray binaries with identified counterparts in which the nature of the outbursts is not clear yet (all these systems are 
discussed in the text). Basic orbital parameters are listed. 

Orbital periods marked with '*' represent the recurrence time of X-ray outbursts; other orbital periods are derived 
from exact orbital solutions. Spectral types have only been included for those systems in which they are derived from 
high signal-to-noise ratio spectra in the classification region. Objects in which the eccentricity is marked as 'large' 
have no orbital solutions, but available data imply e > 0.5. References are only given if data were not included in 
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93.5 


247.5* 


0.66*^^ 


4U 1145-619 


V801 Cen 


B0.2III'' 


292 


187* 


large 


4U 1258-61 


V850 Cen 


BO. 7V'' 


272 


132.5* 


large 


2S 1417-624 


star 




17.6 


42.1 


0.45 


XTE J1946+274 


star 




15.8 


172' 




EXO 2030-f 375 


star 




41.7 


46.03 


0.41^ 


XTE J1543-568 






27.1 


76.6 


< 0.03'' 


2S 1845-024 






94.3 


242.2 


0.88' 


GRO J2058+42 






198 


110*" 




A 0535-668 


star 


B0.5Iir 


0.07 


16.7* 


large 


A 0726-26 


V441 Pup 


08.5V 


103.2 


35* 




A 1118-616 


Wray 977 


09.5V 


406.5 






Cap X-4 


star 


BIV? 


66.3 






SAXJ2239.3+6116 


star 






262*° 





Negueruela & O kazaki ( |200lD 
'' Corbet & Peele fcOOOl) 



^ Negueruela et al. ( |l999| ) 
= Israel et al. (pOOOl) 



° Reig et al. ( |200lD 
* Wilson et al. (|l997|) 



^ M. Scott, pri v. com m. 
j Wilson et al. (|200l|) 
" Wilson et al. (^998|) 



Negueruela, in p rep. 
in't Zand et al. ( |200l|) 
" Negueruela & Coe {^00l\) 



' Wilson et al., in pr ep. 
' Finger et al. (Il999[) 
°in't Zand et al. ( |200C| ) 



should be taken as an upper limit of the distance over 
which the disc outer radius can spread out during one 
orbital period, because the outflow velocity in Be discs is 
certainly much smaller than a few kms~'^. 

From Fig. 0and Fig. ||(a), we note that the truncation 
radii of the Be disc in V 0332-1-53 for a ^ 0.60 is smaller 
than the critical lobe size at periastron. In other words, 
the Be disc never fills the critical lobe even at periastron 
unless the viscosity is very high {a ^ 0.6). 

Table | gives the gap size Ar = rcrit — i^t, where 
?'crit is the mean radius of the critical lobe at perias- 
tron, together with stellar parameters adopted and the 
resulting truncation radii. Table |^ also gives lower lim- 
its for the drift timescale, (Tdrift/^'orb)min, which is given 

by (Tdrift/-Porb)min ~ Ar/[('(;r)max-Porb] ^ Ar/ (0. ICsForb) • 

For a system with (rdrift/^'orb)min > li the gap is wide 
and the truncation is efficient, while the spread of the 



disc can make the truncation inefficient for a system with 

(Tdrift/^orb) mill < 1. 

Given that we expect the viscosity in the disc to be of 
the order of a 0.1 — 1, the Be disc in V 0332-1-53 is likely 
to be truncated at the 3:1 resonance radius or the 4:1 res- 
onance radius. If the disc is truncated at the 3:1 radius, 
which is close to the mean critical radius at periastron 
[Ar/a ~ 0.03 and (Tdrift/-Porb)min — 0.56], a small pertur- 
bation easily causes the outermost part of the disc to fall 
into the gravitational well of the neutron star. However, 
the X-ray history of this system described above suggests 
that the system shows no Type I X-ray outburst in its nor- 
mal state. Therefore, we expect that the disc in V 0332-1-53 
has a viscosity parameter (slightly) less than 0.1 and/or a 
temperature (slightly) lower than ^Tcs- Under those con- 
ditions, it will be truncated at the 4:1 resonance radius 
[Ar/a ~ 0.11 and (rdrift/-Porb)min ^ 2.1]. 
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• (c) EX02030+375 




-2 -1 



x/ a 



- (d) 2S1 41 7-624 




Fig. 2. Description of the orbital model 
adopted for each system, in the refer- 
ence system centred on the Be star. The 
dash-dotted hne represents the orbit 
of the neutron star. The thick dashed 
line represents the position of the first 
Lagrangian point (Li) around the orbit. 
The solid thin lines are the critical lobes 
of the two stars at apastron and perias- 
tron (the position of the neutron star is 
marked with a cross) . The labelled solid 
lines represent the locations of the n : 1 
commensurabilities of disc and binary 
orbital periods at which truncation oc- 
curs for a = 0.3, 0.1, and 0.03 (starting 
from the outside). The distance scale 
corresponding to O.lCs-Porb is shown at 
the lower-right corner of each panel. 



5.2. A 0535+262 

This pulsar has an orbital period Porb = 110.3 d and an 
eccentricity e = 0.47 (Finger et al. 1996a). There is some 



discussion about the exact spectral type of the optical 
counterpart. Most authors sup port a giant in the range 
09.5 -BO, but Wang & Gies (|l998|) cannot rule out a 
main-se quence classification based on UV spectra. Wang 
& Gies ( 1998 ) constrain the mass of the optical compo- 
nent to the range 8Mq < < 22Mq. If the object is a 
giant, the mass should be close to the upper limit. Taking 
= 20Mq, the mass function indicates sini ~ 0.43, 
resulting in an orbital separation a — 286i?Q. 

Four Type II X-ray outbursts from A0535-f262 have 
been observed (in 1975, 1980, 1989 and 1994). The system 
also displays series of Type I outbursts interspersed with 



periods of quiescence. The length of the series or of the 
quiescent states is very variable. A series of three Type I 
outbursts was observed by BATSE during 1993, followed 
by a Type II outburst and later two weaker Type I out- 
bursts (1996a). After that, the source has been in quies- 
cence and the optical counterpart has lost and reformed 



its circumstellar disc (Haigh et al. 1999). 

Fig. shows that the Be disc in A 0535-f 262 is ex- 
pected to be truncated at the 4:1 resonance radius (rt/a ~ 
0.39) for 0.11 ^ a ^ 0.40 and at the 5:1 resonance radius 
(rt/a ~ 0.33) for 0.038 a ^ 0.11. As shown in Fig. |(b) 
and Table || the 4:1 resonance radius is slightly larger 
than the mean radius of the critical lobe at periastron, 
while the higher resonance radii are significantly smaller 
than the critical lobe. The X-ray history of A 0535-1-262 
suggests that the Be disc in this system has a ^ 0.1 and 
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Table 2. Stellar parameters adopted for the systems modelled (those with known orbital parameters) and resulting 
critical lobe radius rcrit, gap size Ar, and lower limit for T^nu/ Poib for different values of the viscosity parameter. 















a = 


0.03 


a — 


0.1 


a = 


0.3 


Name 


Spectral 


M, 










^ Tdrift 


Ar 1 


' Tdrift \ 


Ar / 


' Tdrift \ 


Type 




(Re) 


(K) 


a 




P j 

J^orb / mill 




P I 

\ -'orb / mill 




P j 

^ -'^orb / mill 


4U 0115+634 


BOV 


16.0 


8.0 


3.0 W 


0.46 


0.13 


2.8 


0.08 


1.7 


-0.003 




V 0332+53 


08.5V 


20.0 


8.8 


3.4 10" 


0.50 


0.11 


2.1 


0.03 


0.56 


0.03 


0.56 


A 0535+262 


BOIII 


20.0 


15.0 


3.1 lO" 


0.38 


0.09 


1.2 


0.05 


0.66 


-0.004 




EXO 2030+375 


09III 


22.0 


15.0 


3.4 10" 


0.43 


0.09 


1.7 


0.04 


0.72 


-0.04 






BOIII 


20.0 


14.0 


2.8 10" 


0.42 


0.09 


1.7 


0.04 


0.69 


0.04 


0.69 




BOV 


16.0 


8.0 


3.0 10" 


0.41 


0.12 


2.0 


0.08 


1.4 


0.03 


0.45 


2S 1417-624 


09V 


20.0 


9.0 


3.4 10" 


0.40 


0.10 


1.8 


0.06 


1.1 


0.01 


0.18 




BIV 


12.0 


7.0 


2.5 10" 


0.37 


0.11 


1.9 


0.08 


1.4 


0.04 


0.75 


GRO J1008-57 


BOIII 


20.0 


14.0 


2.8 10" 


0.25 


0.04 


0.38 


0.002 


0.020 


-0.02 






BOV 


16.0 


8.0 


3.0 10" 


0.24 


0.04 


0.41 


-0.004 




-0.03 




2S 1845-024 


BOV 


16.0 


8.0 


3.0 10" 


0.084 


0.014 


0.13 


0.001 


0.013 


-0.010 






BIV 


12.0 


7.0 


2.5 10" 


0.081 


0.015 


0.15 


0.004 


0.038 


-0.007 





^ Td = ^Tcff is adopted for all models. 

^ Vr = O.lCs is adopted as the upper limit on the radial velocity. The value is given for Ar > 0. 



Td ^ and is truncated at the 4:1 resonance radius, 

which will allow enough mass accretion on to the neutron 
star at every periastron passage to cause a Type I out- 
burst regularly (see scenario A in Fig. ^). If this is the 
case, a slight decrease in viscosity (and/or disc tempera- 
ture) will lower the disc radius to the 5:1 resonance radius 
and therefore change the system to a dormant state, as 
has been repeatedly observed. 



5.3. EXO 2030+375 

This transient X-ray pulsar was discovered during a Type 
II outburst in 1985. Since then, it has displayed tens of 
Type I X-ray outbursts in long series separated by periods 
of quiescence. The orbit is characterised by Porb = 46.0 d 
and eccentricity e — 0.41 (Wilson et al. 2001). The opti- 



cal component of this system is heavily obscured and no 
determination of the spectral type has been possible, but 
both the measured mass function and the infrared spec- 
trum imply a spectral type earlier than about Bl (Reig 
et al. 1998). Coe et al. (198S) showed that the colours of 
the object are compatible with a BO spectral type. The 
reddening E(B — V) ~ 3.7 derived by Motch & Janet- 
Pacheco ( 1987) and the distance of 5.3 kpc derived by 
Parmar et al. ( 1989 ) from the change rates in spin period 
and X-ray luminosity can only be compatible if the object 
is a giant. However, given that the distance determination 
has a large uncertainty, the possibility that the object is 
a main-sequence star cannot be ruled out. 

If the spectral type is BOIII, a mass of around M* = 
23M0 is expected. We have adopted a model with M* = 
20Mq and i?* = 14i?f7i, in which case the mass func- 



tion implies an inclination angle i = 56° - this is consis- 
tent with the fact that t he sh ape of the Ha emission lines 
shown by Norton et al. ( 1994 ) is typical of a Be star with 
moderate inclination. The orbital separation would then 



be a = 150i?Q. In order to investigate the dependence of 
our model on the mass of the primary, we have also con- 
sidered the case of — 16 Mq, i?* = 8i?Q (a typical 
BOV star). 

Fig. |l| shows that the difference in ofcrit between the 
BOIII star model and the BOV star model is small. In the 
BOIII star model, the disc is truncated at the 4:1 resonance 
radius {n/a ~ 0.39) for 0.098 <, a <> 0.37 and at the 5:1 
resonance radius [rt/a ~ 0.33) for 0.030 ^ a ^ 0.098, 
while in the BOV star model the disc is truncated at the 
4:1 resonance radius for 0.12 ^ a ^ 0.44 and at the 5:1 
resonance radius for 0.036 ^ a ^ 0.12. We have also tried 
a slightly earlier spectral type, 09III, for the primary, ex- 
pecting that a higher disc temperature would result in a 
larger truncation radius. However, the resulting radii are 
almost the same as those for a BOIII primary, unless a is 
as high as 0.3 (see Table ||). 

As shown in Fig. ||(c), the 4:1 resonance radius is close 
to the mean radius of the critical lobe at periastron, while 
the higher resonance radii are significantly smaller than 
the critical lobe. Our result suggests that the observed 
X-ray behaviour of EXO 2030+375, which has regularly 
exhibited Type I X-ray outbursts, favours a BOIII or a 
09III primary star with a disc with a ^ 0.1 and/or Td ^ 
^Teff that is truncated at the 4:1 resonance radius (see 
scenario A in Fig. 

5.4. 2S 1417-624 

28 1417—624 was observed once in 1978 and was not de- 
tected again until 1994, when a Type II outburst was ob- 
served by BATSE (Finger et al. |l996bD . This was followed 
by five Type I outbursts peaking near apastron. 

The pulsar in this Be/X-ray transient has Pprb — 42.1 d 
and an eccentricity e = 0.45 (Finger et al. 1996b). The 



optical component of the system was studied by Grindlay 
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A. Overflow Induced by Viscous Decretion 

Regular Type I Outbursts 
(a) (b) (c) 




B. Overflow Induced by a Density Wave 
Temporary Type I Outbursts 
(0) (b) (c) 




Fig. 3. Scenarios for two families of Type I outbursts. 
A Regular Type I outbusts induced by viscous decretion 
in an axisymmetric disc, (a) Near periastron, the gas con- 
tinues to overflow through the Li point, replenishing the 
accretion disc around the neutron star. [The accretion disc 
begins to form at panel (d).] (b) At the phase at which 
the critical lobe radius becomes larger than the disc ra- 
dius, the overflow stops and the accretion disc begins to 
fade, (c) While the neutron star orbits far from periastron, 
the disc material drifts outward by viscous decretion. (d) 
At the phase at which the critical lobe radius becomes 
smaller than the disc radius, the gas in the outermost part 
of the disc begins to overflow. The accretion disc is formed 
around the neutron star. B Temporary Type I outbursts 
induced by a slowly precessing density wave, (a) If the 
disc is elongated roughly toward the periastron by chance, 
the overflow occurs around periastron. (b) At the phase 
at which the critical lobe radius becomes larger than the 
elongated disk size, the overflow stops and the accretion 
disc begins to fade, (c) Same as panel (c) in scenario A. 
(d) If the axis of elongation is still roughly toward the 
periastron, at the phase at which the critical lobe radius 
becomes smaller than the elongated disk size, the gas in 
the outermost part of the disc begins to overflow. The 
accretion disc is formed around the neutron star. 



et al. (1984). Owing to the low signal-to-noise of their 
spectra, they did not assign an spectral type, though sev- 
eral He II lines seem to be present in the spectrum dis- 
played, suggesting an 0-type star. Therefore we will con- 
sider again two models, one corresponding to an 09V star 
(M* = 20 Mq, i?* — 9Mq) and one roughly correspond- 
ing to a BIV star (M* = 12 Mq, = 7 Mq). 

As shown in Fig. |l|, the 09V star model and the BIV 
star model give rather different values of Ofcrit: the 09V 
star model results in a value for Ofcrit similar to that in 
the model for A 0535-1-262, whereas the BIV star model 
gives a significantly higher value of acrit- In the 09V star 
model, the disc is truncated at the 4:1 resonance radius 
(rt/a ~ 0.39) for 0.14 ^ a ^ 0.52 and at the 5:1 resonance 
radius (rt/a ~ 0.33) for 0.047 <, a <, 0.14, while in the 
BIV star model the disc is truncated at the 4:1 resonance 
radius for 0.34 ^ a ^ 1.2 and at the 5:1 resonance radius 
for 0.12 ^ a ^ 0.34. Therefore, if the optical counter- 



part of 2S 1417-624 is a BIV star, the disc is likely to be 
truncated at the 5:1 resonance radius or the 6:1 resonance 
radius, which is significantly smaller than the critical lobe 
at periastron (see also Table |^) . In this case, the system 
is unlikely to display Type I X-ray outbursts in a normal 
state, which is consistent with the X-ray history of this 
system. The 09V star model is also consistent with the 
X-ray history if a ^ 0.1. 

We interpret the sequence of Type I outbursts which 
followed the 1994 Type II outburst as a temporary phe- 
nomenon caused by a strong disturbance in the disc. After 
Type II outbursts, the disc is expected to be strongly dis- 
turbed (see Negueruela et al. 2001, from now on Paper 
II, for a discussion of the possible causes of Type II out- 
bursts). Such a disc is likely to be strongly asymmetric 
and to be temporarily capable of fuelling the neutron star 
(see scenario B in Fig. ||). 

5.5. GROJ1008- 57 

BATSE discovered GRO J1008-57 during a Type II out- 
burst in July 1993. Four weak X-ray outbursts separated 
by « 248 days followed. The source is believed to emit 
some X-ray flux during quiescence, but no further out- 
bursts have been observed. The optical counterpart was 
identifled by Coe et al. ( 1994 ), but no information exists 
about its exact spectral type. 

The orbit of this pulsar is not exactly determined. 
However, from the analysis of data from the BATSE ex- 
periment, a most probable orbit can be determined if the 
separation between Type I outbursts (247.5 d) is taken to 
be the orbital period (M. Scott, priv. comm.). In that case, 
the eccentricity is e = 0.66 and Ox sinz = 668.0 ± 9 It — s. 
Again the colours of the primary are compatible with a BO 
star. Therefore we use models for a BOV (M* = 16 Mq, 
R* = 8Rq) and a BOIH (M* = 20 Mq, R^ = URq) 
primary. 

Fig. 1^ shows that both stellar models give similar val- 
ues of acrit- In the BOIII star model, the disc is truncated 
at the 7:1 resonance radius {rt/a ~ 0.27) for 0.16 ^ a ^ 
0.36 and at the 8:1 resonance radius {n/a ~ 0.24) for for 
0.082 ^ a ^ 0.16. In the BOV star model, acrit becomes 
higher by about 20%. 

From Fig. ^ (e), we immediately observe that the trun- 
cation radius is close to or slightly beyond the critical 
lobe radius at periastron unless the viscosity is very low 
(a ^ 0.03). It should be noted that this is a typical fea- 
ture in systems with high orbital eccentricity, because the 
resonance radii are distributed more densely for higher res- 
onances. In other words, disc truncation is less efficient in 
systems with high orbital eccentricity than in systems with 
low or mild eccentricity. Moreover, GRO J 1008— 57 has a 
rather long orbital period. Consequently, the Be disc in 
this system can spread out significantly while the neutron 
star is far away from the Be star. These features should 
enable GRO J1008— 57 to display Type I outbursts regu- 
larly. 
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5.6. 25 1845- 024 

2S 1845—024 has been detected by several satellites at lu- 
minosities compatible with Type I outbursts (see Finger 
et al. 1999 for references). Between the launch of the 



ComptonGRO satellite and 1997, BATSE detected Type I 
outbursts from 2S 1845—024 at every periastron passage. 
This behaviour has continued at least until 2000 (Finger, 
priv. comm.). 

This X-ray pulsar is in a very eccentric (e = 0.88) and 
wide (Porb = 242.2) orbit (Finger et al. |1999D . The optical 
counterpart has not been identified, but the mass function 
constrains it to have Af* ^ 7 Mq. Given that it has an ex- 
act orbital solution, we have modelled the system in spite 
of the lack of data on the optical companion. We have 
considered two models: one in which the primary is our 
typical BOV star (A/, — 16A/q, i?, = 8 Rq) - adopting 
the parameters of a BOIII star did not alter significantly 
the value of acrit - and one with — 12.0 M^, which is 



close to the value used by Finger et al. (|1999| ) for their evo- 
lutionary model (Af* = 11.3 ATq), and then i?,* = 7 Rq 
for a typical BIV star. 

As in GRO J1008-57, both steUar models make Httle 
difference in the value of ckcrit (Fig. |l|), and the truncation 
radius is close to or slightly beyond the critical lobe radius 
at periastron unless the viscosity is very low (a ^ 0.03) 
[Fig. H(f) and Table ^. With the extremely high orbital 
eccentricity, the disc truncation in this system is expected 
to be even less efficient than in GROJ1008— 57. In addi- 
tion, the system has a long orbital period. Therefore, it is 
no surprise that 2S 1845—024 has regularly shown Type I 
outbursts without failure. 



6. Discussion 

In Paper I, we analysed the disc size of V635 Gas, the op- 
tical counterpart of 4U 0115-1-63, and found that the trun- 
cated disc size is significantly smaller than the distance 
to the first Lagrangian (Li) point at periastron. Among 
our conclusions, we found that it is precisely this wide 
gap between the disc outer radius and the position of the 
Li point that prevents the system from displaying Type I 
outbursts, by drastically reducing the accretion rate on 
to the neutron star. In Paper II we pointed out that, for 
the very same reason, the density in the disc of V635 Gas 
is likely to grow with time, to a point where the disc be- 
comes optically thick and unstable to the radiation-driven 
warping, which seems to be at the origin of Type II X-ray 
outbursts. 

As discussed in the end of Sect. ^, we cannot expect 
a perfect truncation. Since the disc density is expected 
to decrease rapidly beyond the truncation radius, the gap 
size determines how effective the disc truncation is. For 
systems with a wide gap, like 4U 0115-1-63, the truncation 
is so effective that the system does not show Type I out- 
bursts under normal conditions. Whatever little amount 
of disc material manages to overcome the tidal trunca- 
tion may be easily ejected from the vicinity of the neu- 



tron star by the propeller effect (Stella et al. 1986). On 
the other hand, for systems with a narrow gap, the trun- 
cation will not be efficient, allowing the neutron star to 
capture gas from the disc at every periastron passage and 
display Type I outbursts. 

Therefore, the main purpose of this paper is to find 
out the general trend in the dependence of the gap size 
on various parameters. A second purpose is to understand 
why we can see two different kinds of Type I outburst 
series: regular series of outbursts of similar intensity in 
some systems and short series of outbursts of varying in- 
tensity (generally associated with Type II outbursts) in 
other systems. 

6.1. Systems with large eccentricity 

We have found that, in systems with high orbital eccen- 
tricity, say e ^ 0.6, such as GRO J1008-57 (e = 0.66) and 
2S 1417-624 (e = 0.88), the truncation radius of the Be 
disc is very close to, in the sense that (Tdrift)min ^ -Porb) 
or slightly beyond the critical lobe radius at periastron 
unless the viscosity is very low (a ^ 0.03). Under such 
conditions, we expect disc truncation not to be efficient, 
allowing the accretion by the neutron star of enough mass 
at every periastron passage to cause a Type I X-ray out- 
burst regularly (see scenario A in Fig. ^). Moreover, the 
long orbital period of GROJ1008-57 (Porb = 247.5 d) 
and 2S 1845-024 (Po^b = 242.2 d) will enable the Be disc 
in these systems to spread out significantly while the neu- 
tron star is far from periastron. This larger drift time will 
allow the neutron star in such systems to accrete rather 
more matter at periastron passage than in a system with 
shorter orbital period. 

Among systems which we expect to be relatively well 
described by these conditions, we list 4U 1258—61 (GX 
304—1) which displayed regular Type I outbursts every 
132.5 d during most of the 1970s (Corbet et al. |1986| ) 
and which has been dormant ever since because of the 
disappearance of the Be disc. Similarly, the behaviour of 
4U 1145— 619, which displays short and not very strong 
outbursts every « 188 d, is indicative of a large eccentric 
orbit. This is supported by the lower limit e ^ 0.6 found 



by Gook & Warwick (1987) from the analysis of changes in 
the X-ray pulse timing. Occasionally, some of the periodic 
outbursts from 4U 1145—619 are rather stronger (challeng- 
ing the conventional definition of Type II outbursts, since 
they cannot be described as very bright). This may be 
related to a large global perturbation in the Be di sc, re - 
flected in the profile of emission lines (Stevens et al. 1997). 
RX J0812.4— 3114 has also shown a long series of Type I 
outbursts without any intervening Type II outburst dur- 
ing a 4- year qua si-cyc le of Be disc formation and dissipa- 
tion (Reig et al. |200lD . Finally SAX J2239.3-h6116, which 
has shown several Type I outbursts separated by w 262 d, 
could also belong to this category. 

One very peculiar system is A 0535—668, which has a 
very short orbital period (Porb = 16.7 d) and is believed to 
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have a very high eccentricity (Charles et al. 1995). If the 
eccentricity is really e 0.8, the neutron star passes close 
to the surface of the Be star at periastron. Our model is 
simply not applicable to such an extreme system. 

It must be noted that the observed behaviour of 
GRO J1008— 57 is only partially compatible with the or- 
bital model. The Type I outbursts are short and not very 
strong, with similar luminosities, as expected for such a 
high-eccentricity system. However, the fact that only four 
have been observed and that they occurred after a Type II 
outburst casts some doubts about the validity of the or- 
bital solution used. 



6.2. Systems with low eccentricity 

It is well known that accretion discs in circular binaries 
with mass ratio 0.05 ^ q ^ 0.25 are truncated at the 3:1 
resonance radius by the tidally driven eccentric instability 
(e.g., Osaki |1996| ). Since the mass ratios of all Be/X-ray 
binaries fall within this range and the mechanism of trun- 
cation is the same for decretion discs as for accretion discs, 
we expect that Be discs in systems with very low orbital 
eccentricity, say e ^ 0.2, are also truncated at the 3:1 
resonance radius. Consequently, the gap size in these sys- 
tems should be much wider than in systems with larger 
eccentricities. Wide gaps in low-eccentricity systems will 
result in disc truncation being so effective that no Type I 
outbursts should occur unless the Be disc is very strongly 
disturbed. Therefore it is likely that systems with very low 
orbital eccentricity will show only Type II outbursts (and 
perhaps temporarily Type I outbursts only when the disc 
is strongly disturbed). 

From this point of view, GS 0834—430 is an extraor- 
dinary system, since it has rather low eccentricity (con- 
strained to be 0.10 ^ e ^ 0.17), but has shown a long 
series of Type I outbursts (Wilson et al. 1997). It must be 
noted that the abrupt change in the separation between 
outbursts observed during this series is highly unusual. 
Moreover, the decrease in peak intensity of the outbursts 
along the series suggests that the Be disc was strongly 
perturbed at the time of the outbursts. We have to admit, 
however, that with the information available, our model 
obviously fails to explain this system. 

Apart from X Per, which shows no outbursts, sys- 
tems with low eccentricity are XTE J1543-568 (e < 0.03), 
which has been detected only recently and does not seem 
to show any clear modulation of its X-ray lightcurve at 
the orbital period (in't Zand et al. |200l| ), and 2S 1553-542 
(e < 0.09), which was only detected once in 1975 (Kelley 
et al. 1983). Our interpretation of their behaviour is that 
these objects, because of effective disc truncation, never 
display Type I outbursts. Their discs therefore accumulate 
mass until very large perturbations develop. Only occa- 
sionally will one of these perturbations result in transfer 
of matter towards the neutron star, producing a bright 
Type II outburst. 



In this respect, it is necessary to insist that the physi- 
cal mechanisms leading to Type II outbursts are not well 
understood. Our model indicates that the truncated discs 
in Bc/X-ray binaries cannot reach a steady state and will 
continuously grow denser. As a consequence, we expect 
the discs to become dynamically unstable. In the case of 
4U 0115-1-63, which has been carefully studied (and per- 
haps also in A 0535-)-262), there is observational evidence 
that links big dynamical instabilities and Type II out- 
bursts (Paper II). From this observational fact, we con- 
clude that the dynamical perturbation somehow results in 
the transfer of large amounts of material from the Be disc 
to the neutron star. No aspect of our model implies that 
such mass transfer must occur. It may well be that in sys- 
tems with low eccentricity, dynamical instability generally 
leads to the collapse of the disc and its fallback on to the 
Be star without any substantial amount of material reach- 
ing the neutron star. If this is the case, low-eccentricity 
transients would be completely undetectable most of the 
time. 

Until recently, it had been generally assumed that the 
orbits of Be/X-ray binaries were in general rather eccentric 
owing to the effect of supernova kicks. The recent deter- 
mination of orbital parameters for several systems with 
low eccentricity (Delgado-Marti et al. 2001; in't Zand et 
al. 2001) may represent a challenge to this idea. If there is 
a selection effect against the detection of systems with low 
eccentricities, their population could be rather larger than 
previously suspected. If eccentric systems display X-ray 
outbursts rather more frequently than systems with low 
eccentricities, they would be more likely to be detected 
over a limited time-span (such as the time since the start 
of X-ray astronomy). 

The orbit of the pulsar GRO J1948-I-32 is constrained 
to have an eccentricity e ^ 0.25 and an orb ital period 
35 d < Porb < 70 d (Chakrabarty et al. |l995| ). Its recent 
identification with the Be/X-ray transient KS 1947-1-300 
means that it has probably displayed at least three Type 
II outbursts between 1989 and 2001. An accurate deter- 
mination of its orbit would therefore be very important 
in understanding whether this selection effect really ex- 
ists. If its eccentricity turns out to be very low, then we 
will have to conclude that low-eccentricity systems may 
display Type II outbursts as frequently as systems with 
moderate eccentricities. If, in contrast, its eccentricity is 
close to e = 0.25, the available data would suggest that 
systems with low eccentricities have less frequent Type II 
outbursts. 

A second system that could serve to understand the 
outbursting behaviour of systems with low eccentricities 
is the LMC transient EXO 0531.1-6609. Based on a sta- 
tistical parameter study, Dennerl et al. (1996) conclude 
that this pulsar (Pg = 13.7 s) is likely to have an or- 
bital period Porb — 25.4 d and low eccentricity e « 0.1. 
EXO 0531.1— 6609 has displayed at least three relatively 
bright outbursts extending over more than one orbital pe- 
riod (i.e.. Type II) in 1983, 1985 and 1993, but also seems 
to have been occasionally detected by ROS AT at a low X- 
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ray luminosity close to periastron. Obviously confirmation 
of the orbital parameters is necessary before any conclu- 
sions can be drawn. 



6.3. Systems with moderate eccentricity 

Four systems (V 0332-1-53, A 0535-1-262, EXO 2030-^375, 
2S 1417— 624) among the six Be/X-ray binaries discussed 
in the previous section have moderate eccentricity. As we 
have already seen, the model results for these systems turn 
out to depend on rather subtle details, in contrast to the 
much more robust predictions for systems with high ec- 
centricity and very low eccentricity. Since the spacing of 
the candidate truncation radii is rather larger in these sys- 
tems than in high-cccentricity Be/X-ray binaries, a small 
difference in the system parameters (orbital, stellar or disc 
parameters) can produce significantly different X-ray be- 
haviour, depending on whether the resulting truncation 
radius happens to be close to the critical lobe radius or 
not. Therefore we have not tried to present general pre- 
dictions for these systems. Instead, in the previous sec- 
tion we have estimated the truncation radii and system 
parameters by comparing the model results with the ob- 
served X-ray behaviour. In this sense, the current status 
of our model can still be considered to be "at the qualita- 
tive level" for mildly eccentric systems. We estimated that 
Be discs in V 0332-^53, A 0535-1-262 in the X-ray active 
state, and EXO 2030-1-375 are truncated at the 4:1 reso- 
nance radius, while those in A 0535-1-262 in quiescence and 
2S 1417— 624 are truncated at the 5:1 resonance radius. 
The smaller truncation radius for A 0535-1-262 in quies- 
cence might suggest a slightly lower disc temperature in 
this state than in the X-ray active state. To understand 
more about these systems, we need a more sophisticated 
model for the interaction between the Be disc and the 
neutron star. 

It is worth noting, however, that all the intermediate 
eccentricity Bc/X-ray binaries known have displayed Type 
II outbursts. This is the case for the four systems listed 
here and also for GRO J1750— 27, a presumably very dis- 
tant source (with Porb = 29.8 d and e = 0.36) which has 
only been detected once during a Type II outburst (Scott 
et al. 1997| ) . This fact suggests that disc truncation is gen- 
erally rather effective for orbital eccentricities e ^ 0.5, 
even though the discs in A 0535-f 262 and EXO 2030-^375, 
under certain conditions, seem to reach a state in which 
truncation is not so effective and long series of Type I out- 
bursts occur (not surprisingly, these two systems have the 
largest eccentricities and longest periods among the five 
known) . 

Though systems with moderate eccentricities domi- 
nate the sample of Be/X-ray binaries for which orbital 
solutions exist at present, it is important to note that at 
least two important selection effects may have been at 
work in the definition of this sample: 



systems with high eccentricities and perhaps also more 
likely than systems with low eccentricities (see above) . 
Their detection as new X-ray sources is therefore eas- 
ier. 

Among Be/X-ray binaries, 4U 0115+63 has been 
by far the most active during the era of X-ray astron- 
omy, displaying no fewer than 13 Type II outbursts 
between 1969 and 2000. From the estimates in Paper 
II, we come to the conclusion that the enhanced activ- 
ity of 4U 0115-1-63 is not owing to causes intrinsic to 
the Be star (such as a particularly high mass-loss rate) 
and therefore it is likely to be connected to its orbital 
parameters. 

The combination of moderate eccentricities and rel- 
atively close orbits results in rather efficient disc trun- 
cation, which allows the storage of matter in the Be 
disc, but at the same time results in gaps which are 
not too wide and may permit the transfer of part of 
this material to the neutron star under certain condi- 
tions. 

2. Once the system has been detected, a relatively close 
orbit may be solved from X-ray observations spanning 
one single Type II outbursts, while in a larger orbit 
a much longer baseline will be needed (only recently, 
permanent coverage over long time-spans with BATSE 
and RXTE has allowed the solution of some orbits with 

Porb ^ 50 d). 

Because of these two factors, judgement on what sort 
of X-ray behaviour is more "typical" among Be/X-ray bi- 
naries has to be suspended until a rather larger sample is 
known. 



6.4. Predictions 

Many of the systems under discussion have shown short 
series of irregular Type I outbursts. This category in- 
cludes 4U 0115+63, V 0332+53, 2S 1417-624 and presum- 
ably GRO J2058+42, for which no orbital solution exists. 
These series have generally occurred in connection with 
Type II outbursts (in most cases, after a Type II out- 
burst). The behaviour of A 0535+262 has also been similar 

we expect 



on some occasions. As mentioned in Sect. 5.4 



1. Be/X-ray binaries with moderate eccentricities are 
much more likely to display Type II outbursts than 



that the temporary Type I outbursts occur in systems in 
which the Be disc has a radius significantly smaller than 
(but not too small compared to) the critical lobe radius 
under normal conditions. These systems will be mostly 
mildly eccentric systems. If the Be disc in such a sys- 
tem is strongly disturbed by, e.g., radiation-driven warp- 
ing and/or a global density wave, it will become strongly 
asymmetric. Then, if the disc elongation is roughly in the 
periastron direction, the disc will be capable of fuelling 
the neutron star and cause a sequence of Type I X-ray 
outbursts (see scenario B in Fig. |). Since GRO J2058+42 
has apparently displayed similar behaviour, we expect it 
to have moderate eccentricity. 

A relationship between the X-ray behaviour and the 
orbital eccentricity stems naturally from the truncated 



Okazaki & Negueruela: X-ray outbursts from Be/X-ray binaries 



13 



disc model. In general, the presence of regular Type I 
X-ray outbursts indicates that the truncated disc around 
the Be star has a radius close to or slightly larger than 
the critical lobe radius at periastron. Such a situation oc- 
curs in all systems with highly eccentric orbit, some mildly 
eccentric systems, and no systems with low eccentricity 
(GS 0834— 430 may be an exception). In a rough sense, 
the lower the orbital eccentricity, the wider the gap be- 
tween the truncation radius and the critical lobe radius 
and therefore the more difficult the accretion on to the 
neutron star. 

A source which hardly fits in this picture is A 0726—26. 
This persistent low-luminosity X-ray source which has 
never shown outbursts displays, however, a 35-d modu- 
lation in its X-ray emission (Corbet & Peele 1997 ). If this 
is the orbital period, A 0726—26 presents very few similar- 
ities to other Be/X-ray binaries, falls very far away from 
the Ps/PoTh correlation and cannot be explained by our 
model. 

Among the Be/X-ray binaries with known counter- 
parts, two systems have been very rarely detected as X- 
ray sources, and even then they have been relatively weak. 
A 1118—616 has only been detected twice, in 1975 and in 
1991. In both cases, the outburst was long, but not very 
bright. Cep X-4 has been detected four times between 1972 
and 1997. All the outbursts have also been relatively weak. 
The two sources display permanent low-luminosity X-ray 
emission when not in outburst. In both cases, the counter- 
parts present evidence for very large circumstellar discs. 
Therefore we expect these two systems to have wide and 
not very eccentric orbits. 

We note that even at the qualitative level our model 
has some predictive power and can thus be tested. Finding 
out, for example, that RX J0812.4— 3114 has an orbital 
eccentricity close to zero or that A1118— 616 has a very 
eccentric orbit would certainly force us to reconsider most 
of the conclusions reached here. 

7. Concluding Remarks 

We have applied the resonantly truncated disc model de- 
veloped by Negueruela & Okazaki (2001) to six systems 
which have displayed Type I X-ray outbursts. We have 
found that the model naturally explains the X-ray be- 
haviour of these systems and that all systems are consis- 
tent with having a similar viscosity parameter a ^ 0.1. 

In our model, regular Type I X-ray outbursts occur in 
all systems with high orbital eccentricity, some mildly ec- 
centric systems, and no systems with low eccentricity. On 
the other hand, systems which have temporarily shown 
Type I outbursts will be mostly mildly eccentric systems, 
in which the Be disc temporarily becomes strongly asym- 
metric. 

In the scenario presented here for Be/X-ray transients, 
disc material can reach the neutron star only via the Li 
point, and will therefore have a not very high velocity rel- 
ative to the neutron star. Since such a flow (which may 
be said to represent the disc version of Roche-lobe over- 



flow) carries high angular momentum, an accretion disc 
may temporarily be formed around the neutron star dur- 
ing each X-ray outburst. Such temporary accretion discs 
have probably been observed in 2S 1845—024 (Finger et al. 



1999), but extensive searches have failed to detect them 



around A 0535-1-262 (Motch et al. |199lD . The accretion 
flow in the vicinity of the neutron star (i.e., once it has 
left the effective Roche lobe of the Be star) should be mod- 
elled in detail in order to understand if such discs should 
always form and to predict the X-ray lightcurves derived 
from our scenario. 
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